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A zero temperature Anderson-Mott transition driven by spin disorder can be 'tuned' by an applied 
magnetic field to achieve colossal magnetoconductance. Usually this is not possible since spin 
disorder by itself cannot localise a high density electron system. However, the presence of strong 
structural disorder can realise this situation, self consistently generating a disordered magnetic 
ground state. We explore such a model, constructed to understand amorphous GdSi, and highlight 
the emergence of a spin glass phase, Anderson-Mott signatures in transport and tunneling spectra, 
and unusual magneto-optical conductivity. We solve a disordered strong coupling fermion-spin- 
lattice problem essentially exactly on finite systems, and account for all the qualitative features 
observed in magnetism, transport, and the optical spectra in this system. 



The 'Anderson-Mott' insulator-metal transition (IMT) in 
disordered interacting systems [1,2] and the physics of 
'colossal magnetoresistance' (CMR) in, for example, the 
manganese oxides [3] are topics of great current interest. 
Even though the "ultimate" example of CMR would be a 
magnetic field driven, zero temperature, insulator-metal 
transition, such a scenario had not been realised experi- 
mentally till recently. Most observations of CMR are at 
finite temperature, across a ferromagnet to paramagnet 
transition [3-5]. Canonical 'Anderson-Mott' systems do 
not show large magnetoconductance (MC) and the stan- 
dard CMR systems do not involve localisation physics. 
These two fields of research have evolved independently. 

Experiments [6-10] on amorphous a-GdSi reveal that 
the presence of doped magnetic moments in a strongly 
disordered system can combine features of the standard 
doping driven IMT in amorphous systems [11] with the 
physics of field driven IMT and CMR. The magnetic 
ground state in such a system is a spin glass. There 
are distinct signatures of electron correlation in the con- 
ductivity and tunneling spectra, and huge transfer of op- 
tical spectral weight on application of a magnetic field. 
The experimental observations, discussed below, cannot 
be understood within the standard scenarios [1-5] de- 
veloped for IMT and CMR and require an independent 
and comprehensive framework. Our main achievement in 
this paper is (i) to provide the first understanding of the 
unique properties of this system, and (ii) demonstrate a 
many body technique that allows controlled approxima- 
tions in a strongly disordered interacting system. 

The measurements have been made on a-GdSi and 
simultaneously on the non-magnetic analog a-YSi. 

(i) Both Yj-Sii-a; and Gda-Sii-a; show an IMT [6] as the 
doping, X, is increased beyond a critical value, Xc- The 
critical doping Xc ~ 14% in YSi and Xc ~ 15% in GdSi. 

(ii) For a; ^ Xc in GdSi, a magnetic field, hc{x), can ac- 
tually drive an insulating sample metallic. YSi samples 



show weak positive magnetoresistance. (Hi) The density 
of states (DOS) at the Fermi level, iV(0), in GdSi, probed 
through tunneling conductance measurement [7], grows 
as {h—hcY' across the IMT, while adc increases as (h—hc). 
(iv) The optical conductivity in GdSi shows large transfer 
of spectral weight [8] to low frequency from ui > OA eV 
on application of a field of a few Tesla. Transfer of spec- 
tral weight also occurs on raising temperature, and this is 
seen in both GdSi and YSi. (w) The low field a.c suscep- 
tibility in GdSi reveals [9] that the magnetic degrees of 
freedom freeze into a spin glass state at low temperature. 
The freezing temperature, T/, increases from ~ 1 K at 
X = 0.04 to T/ - 6 K at X 0.20. The Y doped samples 
are diamagnetic. (vi) The 'effective moment' inferred [9] 
from x{T) differs from the expected value for Gd, and the 
high temperature magnetic specific heat [10] per doped 
Gd is almost 50% larger than log{2S -\- 1). 

Observation (i) above is standard in disordered sys- 
tems, {ii) — {iv) would be expected in CMR materials, 
and (w) — {vi) seem to be unique to the combination. 

A 'first principles' model for amorphous GdSi will have 
to consider an underlying 'random' structure in which a 
fraction x of the sites are occupied by Gd atoms and 
(1 — x) by Si atoms. The Gd and Si atoms have different 
orbital structure so a complicated set of inter-orbital, in- 
tersite hopping possibilities need to be considered. We 
try to retain the essential features in the following, sim- 
pler, one band model: 



H= -t ^ c\^Cj„ + l^)^i'y ~ X! ^^-^^ 



- A ^ n^x^ + ^K^xl + Hco 



(1) 



We use a tight binding model with uniform hopping t, 
and an on site potential uniformly distributed between 
±A/2. The sites labelled 'i^' are a fraction x of the lat- 
tice corresponding to the dopant (Y/Gd) locations. The 
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electron-spin coupling is J', are the S = 7/2 Gd spins, 
and (Tifi, = X^a^ c|(jT^''cj/3, where T/j, are the Pauli matri- 
ces. The local displacement variables (bond distor- 
tions [12]) coupled to the electron density via A, and the 
structural stiffness is K. Hcoui would include Hubbard 
and long range Coulomb interactions. 

The width of the impurity level distribution (A) in 
YSi/GdSi has been estimated [8] to be ~ 200 meV, and 
the 'polaron binding energy' g = }? jK ~ 30 meV. The 
existence of lattice polaron effects in amorphous semicon- 
ductors had been argued early on by Anderson [12], and 
has been revived now [8] in the context of doped a-Si. 
We think the diamagnetism in YSi confirms bipolaronic 
effects, but these lattice effects are probably not very 
important in GdSi. There is no simple estimate of the 
'effective hopping amplitude' to be used in a single band 
approximation. However, calculations on the Anderson 
model indicate [13] that we need A/t « 14 to localise 
10% of the electronic states in the band. This suggests a 
rather small effective hopping amplitude ^ 200 K, if A is 
200 meV. The electron-spin coupling J' S (called J' from 
now on), arising out of the — / coupling in Gd, is large. 
It is estimated to be ~ 0.9 eV from photoemission mea- 
surements [14] on Gd, but would be somewhat smaller 
in the effective one band description that wc arc using. 
Although the parameter values have some uncertainty it 
is clear that A, J' > t. We use A/t = 11, J' jt = 4, 
g/t = 0.5; roughly consistent with the experimental esti- 
mates. The electron-phonon and electron-spin coupling 
are operative only at the dopant sites. We distribute the 
'impurities' (Y or Gd in the Si host) into Nx sites with 
the lowest potential in any given realisation {ci}. This 
ensures that at low dopant concentration, the electrons 
are trapped near the dopant sites. Wc measure all en- 
ergies in units of t, and finally assume t ~ 500 K for 
comparing our energy scales with the data. 




FIG. 1. (a). Order parameter for freezing, q{T) (see text), 
within our effective magnetic model. Insets: (6). The inverse 
susceptibility x~^{T) over a large temperature range, (c). 
Freezing temperature Tf{x). Simulation on 10^ lattices. 




FIG. 2. (a): D.C conductivity, at T = 0. Model for: YSi 
(circles), GdSi in the spin glass phase (squares), and GdSi in 
the fully polarised phase (triangles). (6): L — > oo extrapo- 
lation to construct (Jdc (sec text), (c): GdSi: dependence of 
Udc, at a; = 0.10 and T = 0, on magnetisation. 

Our principal results using the Hamiltonian above are 
the following {i) the magnetic ground state in the low 
doping region is a spin glass with T/ having overall scale 
^ /{J' + A) and following the experimental doping 
dependence (Fig.l), {ii) there is a metal- insulator tran- 
sition for both GdSi and YSi with decreasing x, with 
^Gd ~ and a field driven IMT and OMR for GdSi 
(Fig. 2), {in) there is large transfer of spectral weight from 
high to low frequency (Fig. 4) in cr(a;), driven by an ap- 
plied magnetic field in GdSi. We explain our scheme of 
calculation next and then discuss these results in detail. 

Since the Hamiltonian involves A/t, J'/t » 1 and 
(7 ~ C>{t) none of these couplings can be handled pertur- 
batively. To study the properties of this model within a 
controlled approximation we use a finite size combination 
of Monte-Carlo and exact diagonalisation [15] (MC-I-ED). 
This approach exactly handles the strong disorder, but 
treats the spin and lattice variables as 'classical'. Since 
we have 25* ^ 1, the 'classical' spin limit should be a 
reasonable starting point. At strong disorder, the lead- 
ing effect of phonons should also be accessible classically. 

If we ignore Hcoul to start with, H represents non- 
interacting fermions coupled to classical variables Si and 
Xi, in addition to the random potential e^. The ej 
are 'quenched' variables while the spin and lattice de- 
grees of freedom are 'annealed', with the distribution 
P{x,S} = Z-^Tre-f^" where Z = J VSVxTre-^" is 
the full partition function for a specific realisation of {cj}. 

The 'exact' MC+ED allows only small system sizes, 
O(IOO) sites, so the key step is to construct an approxi- 
mate 'effective Hamiltonian' for the lattice and spin vari- 
ables. Once the magnetic and phonon problem arc self 
consistently solved, the T = electron problem can be 
solved in the classical ground state {S,, Xjjo, which itself 
depends on {ci}, finally averaging over disorder. 

Formally the magnetic effective Hamiltonian is 
Heff{S} = -i log f VxTre-/^". The magnetic prob- 
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lem involves J' jt^ 1, a dilute system (the spins occupy 
only a fraction x of sites) , and strong disorder in the elec- 
tron system. There is no perturbative expansion possible 
in J' but the large J' and .x <C 1 allows a simplification. 
In this limit the doped carriers are essentially localised 
at the magnetic sites, with the electron density falling 
off exponentially away from the sites. This generates a 
pairwise antiferromagnetic coupling, for Rij > 1, with 
Jij ^ {t^/J')e-^*i/^(-^'\ with A oc 1/J'. For neighbour- 
ing sites there is a ferromagnetic coupling 0{t). 



N(e) 




FIG. 3. (a). Density of states in the model for GdSi for two 
densities, a; = 0.1 and x = 0.2, for random spins (m = 0) and 
polarised spins (m = 1). The curves are vertically shifted by 
0.04. Fermi energy marked by arrows, {h). DOS in model 
for YSi, at a; = 0.1 (above) and x = 0.2 (below), (c). DOS 
at €f in the model for GdSi, at a; = 0.1, with increasing 
magnetisation. 

In the disordered system, for a given separation Rij, 
the bonds have a distribution. We construct the dis- 
tributions P{J,Rij) and study the model: iJe//{S} = 
X^ij Jij^i-Sj- For a pair of moments located at Rj and 
Rj the Jij is picked from P{J,Rij). This ignores corre- 
lations between bonds in a specific {ci} realisation. We 
simulate the model for different dilutions, compute the 
order parameter for freezing, q(T) = (xN)^^ 1(8^)^1, 
and check that the structure factor has no peaks at any 
wavevector Q. Fig.l shows q{T), alongwith x(T) at 
X = 0.2, and the freezing temperature [16] Tf{x). At 
X = 0.2 our Tf ~ lOK, while experimentally T/ ~ 6K. 
We have done the exact MC-I-ED simulation for 4^ sys- 
tems, with the same J', A and x ^ 0.1 and vcrificxi [17] 
that the system freezes into a spin glass, with Tf within 
10% of our result here. 

Having established the existence of a glassy state for 
the {S,} we will simplify the remaining electron-phonon 
problem by assuming the spins to be frozen in an uncor- 
related random manner. 

The effective Hamiltonian for phonons is Hej f{x} = 
— ^ log J T>STre~^^ . At moderate g and strong disorder 



there would be 'frozen' bond distortions in the ground 
state, and we cannot expand about the Xi = state. 
To incorporate this effect we use the lowest order self- 
consistent expansion, i.e, Heff{x} « ^KJ2i ^i+J2i ^i^i 
with Ui = —Xfii, where fii = (rij), computed in the 
electronic ground state. The minimum of i?e//i the 
lattice distortion in the ground state, corresponds to 
Xi = {X/K)ni. The T = problem now corresponds to 
electrons in the background of structural disorder, {e,}, 
coupled to randomly oriented spins with coupling J', and 
density coupled to a phonon field Xi = {\/K)ni, i.e. 



Tjel 
^eff 



Hkin + ^ ^i^i — J' ^ 



fifV'n^n^ (2) 



We solve this problem through iterative ED. The trans- 
port and spectral properties of YSi correspond to g = 0.5 
and J' = 0, while for GdSi g = 0.5 and J' = 4. The 
ED is done for a sequence of sizes 6 x 6 x L, with 
L = 24, 32, 40, 48. Due to the finite size gaps the d.c. con- 
ductivity cannot be directly computed on finite systems. 
We use the Kubo-Greenwood formula to compute the in- 
tegrated optical spectral weight (Ti„t( Aw) = alu})dw, 
disorder average, and invert to obtain the optical con- 
ductivity (t(w). The extent of averaging varies from 
400 — 100 realisations, decreasing with increasing L. We 
track a(oJref',L) with coref oc L~^, setting tOref = 0.08 
at L = 32, and use adc = limL^oo(^{ujreS , L). Fig. 2(a) 
shows this 'dc' conductivity appropriate to YSi and GdSi. 
Our smaller than the experimental values and we 

have not fine tuned parameters to match the data. The 
(T arc in units of ire^ /TiaQ and the typical values shown in 
Fig.2(a) are ~ 0.01. For ao ~ 2l, cTdc ~ 400 (ricm)-\ 
roughly as in experiments [6]. Fig. 2(6) shows the L de- 
pendence of ij{ujref', L), while Fig. 2(c) shows the 'magne- 
toconductance'. We have checked that the conductivity 
in a 'spin glass' backgound shows the same trend as for 
random spins and the numbers match within ~ 20%. 
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FIG. 4. Variation of optical conductivity in the model for 
GdSi, at T = 0, with degree of magnetisation. Note the log 
scale in frequency. Inset: magneto-optical conductivity. 



3 



To our knowledge there are no standard results on elec- 
tron systems combining strong structural disorder and 
strong coupling to dilute magnetic moments. In the well 
studied opposite limit, A/t S> 1, J' /t <C 1, spin flip scat- 
tering actually weakens [1] Anderson localisation, so large 
J'/t is crucial to enhanced localisation in GdSi. We have 
cross checked the trends in cr^c by computing the aver- 
aged Greens function C'^ (r — r') at large A, with in- 
creasing J'. Using TrG = G'^'^(r— r') as the indicator 
of 'delocalisation', we find that at |r — r'\/ao = 6, TrG 
grows with J' upto J' ~ 0.5 and then falls rapidly [17]. It 
recovers quickly as the spins are polarised by an applied 
field, tracking the change in conductivity. 

The DOS in GdSi, at low x, has a broad minimum at 
ep, Fig. 3(a), since J' pulls down states to lower energy. 
This minimum is not related to the Altshuler-Aronov 
'correlation gap' which would be a sharper feature [2] 
near ep (with effects of Hcoui included). In YSi, the ef- 
fect of phonons on the disordered background shows up 
as a sharp dip [1] in the DOS, Fig. 3(6), since it generates 
a short range attractive interaction with Ueff = —)? /K. 

We mimic the effect of finite magnetisation (m), in 
GdSi, by using a spin distribution with finite {Sz)- Fi- 
nite m leads to significant redistribution of weight in the 
DOS, due to the large J', which should be visible in pho- 
toemission measurements. The conjunction of increased 
mobility, and increased DOS near e^, Fig. 3(c), leads to 
the large changes observed in (t(cj), Fig. 4. The 'outer 
scale' in cr(w) is ~ 5t « 0.25 eV, as in the data [8]. 

There are certain experimental features for which the 
'Mott' aspect is essential. These are principally the y/T 
dependence in addT), the ^/uj correlation gap, the T 
driven spectral weight transfer in a{u>), and the excess 
magnetic Cy. Most of these are generic correlation ef- 
fects, well known in other amorphous systems [1,11], and 
unrelated to the magnetic character. 

Let us rc-emphasizc the imiqucncss of the system we 
study. Disorder, electron-spin coupling and clectron- 
phonon interactions are features common, in some form, 
to a-GdSi, Anderson-Mott systems (NbSi, say) and the 
CMR manganites. The crucial differences are: (i) GdSi 
is a strongly disordered 'dilute' magnetic system, with 
strong electron-spin coupling. These features are essen- 
tial to the spin glass behaviour and the consequent IMT 
and CMR. Electron-phonon interactions, even if present, 
are not crucial to the physics, {ii) Anderson-Mott sys- 
tems are also strongly disordered, but nominally non mag- 
netic. There are no remarkable magnetic field effects and 
the physics is controlled by disorder and electron correla- 
tions. {Hi) Most CMR manganites are reasonable metals 
at low temperature, indicating weak intrinsic disorder. 
They have strong electron-spin coupling on a periodic Mn 
lattice, which, in contrast to 'dilution', promotes dou- 
ble exchange ferromagnetism. Electron-phonon (Jahn- 
Teller) interactions are important in these systems. The 
finite temperature IMT and CMR are related to multi- 



phase coexistence [3] and not an Anderson transition. 
GdSi differs also from the diluted magnetic semiconduc- 
tors in that the spin polaron concept [5,6] is not tenable 
in this high electron density system, due to strongly over- 
lapping wavefunctions. 

In conclusion, this is the first explanation of insulator- 
metal transition, CMR, spin glass freezing and op- 
tical properties of a-GdSi, bridging the gap between 
Anderson-Mott transition and CMR systems. Our re- 
sults are based on an exact finite size calculation, han- 
dling strong disorder and interactions. The effect of 
Coulomb interactions is understood only qualitatively at 
the moment, and their inclusion would be the next step. 
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